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ABSTRACT 

The solar meridional flow is an important ingredient in Babcock-Leighton type 
models of the solar dynamo. Global variations of this flow have been suggested 
to explain the variations in the amplitudes and lengths of the activity cycles. 
Recently, cycle-related variations in the amplitude of the P\ term in the Legendre 
decomposition of the observed meridional flow have been reported. The result is 
often interpreted in terms of an overall variation in the flow amplitude during the 
activity cycle. Using a semi-empirical model based upon the observed distribution 
of magnetic flux on the solar surface, we show that the reported variations of the 
P\ term can be explained by the observed localized inflows into the active region 
belts. No variation of the overall meridional flow amplitude is required. 

Subject headings: 



1. Introduction 



The solar meridional flow is a large-scale plasma motion which transports material 
from the equator towards the poles near the surface. Temporal variations of t his flow have 



been reported in terms of a reduced flow ve locity during solar maxima (e.g., iKomm et al. 
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Haber et al. 



2008 



2002 
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Gonzalez Hernandez et al.l 120101 ) 



2004 



Basu fc Antial 12003 ) or as local inflows towards the activity belts 



Gizonll2004J : iGizon fc Rempelll2008l : ISvanda et al. 



Recently, lHathaway fc Rightmird ( 120101 ) studied the time dependence of the meridional 
flow by considering magnetic features in SOHO/MDI magnetograms since 1996 as tracers of 
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the flow field near the solar surface. The large-scale motion of the magnetic elements is a 
combination of the real large-scale bulk plasma motions advecting the magnetic features and 
a diffusion-type motion caused by the action of the random granular and supergranular flows 
on l arge-scale gradient s in the distribution of the number density of magnetic field elements 
(see IWang et al.ll2009l ). Rather than disentangling the two effects, lHathaway fc Rightmire 
( 120 lOl ) considered the magnetic field elements to be simple tracers of the large scale flow field 
near the solar surface. 

The projection of the resulting flow profiles on the Legendre polynomial P\ revealed 
a significant time variation of the corresponding coefficient, with lower values during the 
activity ma ximum of cycle 23 an d higher values during the preceding and following activ- 
ity minima. iBasu &: Antial (120101 ) used SOHO/MDI velocity maps for the same period and 
studied the evolution of the meridional flow by helioseismic methods, thereby also obtaining 
information about its variation with depth. The decomposition of the near-surface flows in 
Legendre po lynomials yielded a variation o f the P\ co efficient in qual i tative agreement with 
the result of lHathaway fc Rightmird (j2010l ). However, IBasu fc Antial (120 lOh also found that 
the variation of the meridional flow is connected with a flow pattern that migrates equa- 
torward in parallel with the activity belts as outlined by the butterfly diagram of sunspots. 
This flow pattern probably reflects the latitudinal inflows towards the act ivity belts (e.g., 
Gizon fc Rempelll2008l ; ISvanda et al.ll2008l ; iGonzalez Hernandez et al.ll2010l ). This raises the 
question whether the variation of the P\ coefficient could (partly or fully) be understood as 
the superposition of the local inflows onto an undisturbed large-scale meridional circulation. 

The distinction between a general reduction of the flow speed (which should then also 
affect the return flow in the deep convection zone) and a superposed activity-dependent 
local flow near the surface is important both in terms of its effects and in terms of its 
cause(s). Overall meridional flow variations could possibly play a role in modulating the 
amplitude and length of the activity cylces. In particular, the cycle period of an advection- 
domi nated Babcock-LeiKhton type dynamo is sensitive to the strength of the meridional 
flow (JDikpati fc Charbonneaulll999l). because the fl ow essentially acts like a "conveyor belt" 
in transporting the field ( JDikpati &: Gilmanll2006l ). A localized near-surface inflow, on the 
other hand, does not influence the overall speed of the conveyor belt and thus the cycle 
period. A locali zed inflow does, ho wever, affect the amplitude of the polar field and the axial 
dipole moment (I Jiang et al.ll2010l ). with possible attendant effects on the amplitude of the 
following activity cycle. In terms of the cause of the time variations, theoretical models have 
been sugge sted to explain either t ype of time dependence of the meridional flow (for a recent 



review, see 



Brun & Rempelll2009h . 



The purpose of this paper is to consider whether the variations of the P\ coefficient 



-3- 



reported by lHathaway fc Rightmird (120101 ) and iBasu fc Antial (120101 ) provide conclusive ev- 
idence for a solar-cycle modulation of the overall meridional flow. We find that a semi- 
empirical model of the local near-surface inflows towards the activity belts results in a time- 
dependent P\ signal in the Legendre decomposition that largely explains the observational 
results, without requiring a change in the overall meridional flow amplitude. 



2. Modelling the inflow into the active region belts of cycle 23 

We consider a model of the latitudinal inflows that connects their activity-related vari- 
ation with the measured magnetic field at the solar surface during cycle 23. To this end, 
we assume that the spatial and temporal variations of the inflows are directly related to the 
surface magnetic field distribution. This is consistent with the theoretical model of flows 
towards active regions being driven by the cooling associate d with the excess brightness of 
magn etic features in plage and enhanced network regions (jSpruitl 120031 ; iGizon fc Rempel 
20081 ). The strength of the inflow towards each magnetic feature is taken to be proportional 
to the local unsigned vertical magnetic field, \B\. The flow at any point is the superposition 
of the effect of all magnetic features. The individual inflows are radially symmetric, so that 
the strength of the resulting superposed large-scale flow is proportional to the horizontal gra- 
dient of \B\. Since we are only interested in the axisymmetric component of the meridional 
flow perturbation, we take the speed of the latitudinal inflow, v(X, t), to be proportional to 
the latitudinal derivative of the longitudinally averaged magnetic field, (|S|)^(A, t): the net 
meridional flow perturbation at any location reflects the difference between the effects of 
the inflows driven by the magnetic features lying to the north as against those lying to the 
south. The resulting expression reads 



v(X,t) =c 



d(|g|)»(A,t) 



The constant of proportionality, Cq, is calibrated through matching the resulting flow ampli- 
tudes with those of the observed inflows towards the activity belts. The azimuthal average 
of the radial (vertical) field, determined from SOHO/MD1 synoptic maps u is shown in Fig- 
ure [TJ The data have been remapped from a grid equally spaced in sine latitude to one which 
is equally spaced in latitude (with a spacing of 0.165°). 

We removed small-scale fluctuations from the latitudinal derivatives by smoothing with 
a Gaussian filter with a half width at half maximum of 20°, indicated by the curly brackets 



1 From |http : //soi . stanf ord . edu/magnetic/index6 . html 
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Fig. 1. — Time-latitude diagram of the longitudinally averaged unsigned radial magnetic 
field at the solar surface as derived from SOHO/MDI synoptic maps. 
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in Eq. ([I]). The time- latitude diagram of the resulting meridional flow perturbation is shown 
in Figure [2j The flow polewards of ±60 was set to zero to avoid the problems near the poles 
evident in Figure [TJ which result from the poor determination of the magnetic field near the 
poles and the effect of the varying solar B angle. 

The constant of proportionality in Eq. ([1]) was calibrated as cq = 9.2 ms _1 G _1 deg 
by requiring that the amplitude of the inflow should be comparable to that reported by 
Gizon et al.l (120101 ) . The results given there have, by construction, zero cycle-averaged flow 
at all latitudes. For the comparison we therefore first subtracted from our model flow the 
time-average as a function of latitude. We also took a 5- year running temporal mean since 
the observationally inferred flows are smooth on timescales of years (again by construction). 
These two effects reduce the amplitude of the flo w perturbation to v alues between —5.4 ms _1 
and +4 ms _1 , comparable to those reported by iGizon et al.l ( 120101 ). 



The modeled bulk inflow is into regions of higher unsigned flux densities, and hence it 
drives the magnetic flux elements in the direction of the gradient of the unsigned magnetic 
flux. The sense of the flow is thus opposite to that of the diffusive transport (via a random 
walk) of flux elements away from regions where the number density of such elements is 
high. The nature of the model inflow and diffusive transport implies the need to properly 
dis entagle the two wh en using magnetic elements as tracers of the flow as was pointed out 
bv lWanget~ai1(l2009l ). 



3. Analysis 



To compare the model with the results of lHathaway fc Rightmird (120101 ) and lBasu fc Antia 



(J2010l ). we decompose the modeled flow in Legendre polynomials. There are two main effects 
which cause the inflows towards the activity belts to contribute to the P\ coefficient. The 
first is that the calculation of the coefficient more heavily weights the flow in the range 20° to 
50° than near the equator. As illustrated in Figure EJ the weighting thus favors the latitude 
range where the inflow corresponds to a motion towards the equator. As a consequence, the 
P\ coefficient is reduced with respect to the value for the undisturbed meridional flow. The 
second effect is due to the fact that the inflows extend over latitudinal scales of ~ 20°: when 
the activity belts approach the equator, there is overlapping and partial cancellation of the 
oppositely-directed flow perturbations on both sides of the equator. 

The time variation of the coefficient C2 of the P\ ter m in the Legendre decompositi on is 
shown in Figure HI al ong with the corresponding results of lHathaway fc Rightmird (120101 ) and 
Basu fc Antial (120101 ). In all cases, we subtracted the (different) mean values of the coefficient. 
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Fig. 2.— - Time-latitude diagram of the model meridional flow perturbation, v(X,t) = 
co{d(\B\)^(X,t)/d\}. For a better representation in the figure, the velocities have been 
symmetrized, so that positive velocities correspond to poleward flow in both hemispheres. 
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Fig. 3. — Effect of converging flows towards the activity belts on the amplitude of 
the coefficient to P 2 X in a Legendre polynomial decomposition of the meridional veloc- 
ity. The modeled flow perturbation for July 20, 1999 (solid black curve) is weighted by 
P^(90° - A) cos(A)/ /^(SO - A )] 2 cos(A)dA (dashed blue curve) to obtain the integrand for 
the calculation of the P 2 X coefficient (red dash-dotted curve). The integral of this curve over A 
yields the coefficient of P^- The fact that the near-equator latitudes are attenuated leads to 
a net negative value of the integral over latitude, corresponding to an apparent deceleration 
of the meridional flow. The flow and the weighting function have both been multiplied by 
— 1 (thus leaving the product unchanged) in the south hemisphere (A < 0), in order to make 
clear that both hemispheres contribute to the coefficient of P\ in the same way. 



It can be seen that the mod el inflows reproduce the variati on of the P\ coefficient over the 
solar cycle 23 as reported by lHathaway fc Rightmird (120101 ). including the amplitude of the 
va riation. Fluctu a tions on short time scales are not reproduced. The cycle variation found 
by lBasu &: Antial (120101 ) for a depth of 1.4 Mm is qualitatively similar, but the amplitude is 
about a factor 2 higher - reproducing this would require that our calibration factor, cq, be 
correspondingly higher. 

The observational results differ with regard to a possible difference of the meridional flow 
speed during the two activity minima included in the time series. While lHathaway fc Rightmire 
(120101 ) fin d a somewhat higher speed during the recent minimum in comparison to the pre- 



vious one, iBasu fc Antial (120101 ) do not detect a significant difference. In our model, the flow 
perturbation is related to the surface activity, so that we expect an almost undisturbed flow 
during solar minima and thus no significant differences between the minima. 



Conclusion 



The time variation of the Pq coefficient in the Legendre polynomial expansio n of the 
meridional flow observed by lHathaway fc Rightmird (120101 ) and lBasu &: Antial (120101 ) is qual- 
itatively reproduced by a model describing inflows into the activity belts. The model is based 
upon the instantaneous latitudinal distribution of the unsigned magnetic flux at the solar 
surface. Our results explain the fact that the meridional velocity, determined by the coefi- 
cient of P 2 X m the Legendre expansion, is smaller during activity maxima - without the need 
to invoke a change of the overall meridional flow speed. 

We have modeled the flow for cycle 23 using the observed magnetic fields, which is 
consistent with the idea that the inflow is driven by therm al changes du e to the enhanced 
radiative cooling associated with the plage magnetic fields (ISpruitll2003l ). Both our simple 
model and the underlying theory indicate that the inflows should depend on the strength of 
the activity cycle, so that the ampli tude of the inflows in previous cycles could, in principle, 
be inferred. As recently shown by iJiang et al.l (120 lOl ). such inflows affect the Sun's polar 
fields, with potential implications for dynamo models. 
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Fig. 4. — Time variation of the coefficient of P\ in the Legendre decomposition of the 
meridional velocity for cycle 23. The mean value of the coefficient has been subtracted in 
all cases. The result from the modeled lat itudinal inflows towards the a ctivity belts (red 
curve) is shown toge ther with the results of lfiathawav fc Rightmird ( 2010l black curve) and 
Basu fc Antial (j2010l for 1.4 Mm depth, blue curve). iBasu fc Antial (120101 ) used a different 
normalization, so their result was multiplied by a factor 1.5 to bring it on a common scale. 
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